Selective translational control of cellular and viral mRNAs by RPS3 mRNA binding.

RPS3, a universal core component of the 40S ribosomal subunit, interacts with mRNA at the entry channel. Whether RPS3 mRNA-binding contributes to specific mRNA translation and ribosome specialization in mammalian cells is unknown. Here we mutated RPS3 mRNA-contacting residues R116, R146 and K148 and report their impact on cellular and viral translation. R116D weakened cap-proximal initiation and promoted leaky scanning, while R146D had the opposite effect. Additionally, R146D and K148D displayed contrasting effects on start-codon fidelity. Translatome analysis uncovered common differentially translated genes of which the downregulated set bears long 5'UTR and weak AUG context, suggesting a stabilizing role during scanning and AUG selection. We identified an RPS3-dependent regulatory sequence (RPS3RS) in the sub-genomic 5'UTR of SARS-CoV-2 consisting of a CUG initiation codon and a downstream element that is also the viral transcription regulatory sequence (TRS). Furthermore, RPS3 mRNA-binding residues are essential for SARS-CoV-2 NSP1-mediated inhibition of host translation and for its ribosomal binding. Intriguingly, NSP1-induced mRNA degradation was also reduced in R116D cells, indicating that mRNA decay occurs in the ribosome context. Thus, RPS3 mRNA-binding residues have multiple translation regulatory functions and are exploited by SARS-CoV-2 in various ways to influence host and viral mRNA translation and stability.

Intragenomic rearrangements involving 5'-untranslated region segments in SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses.

BACKGROUND: Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses that was first documented in an infectious isolate of the betacoronavirus SARS-CoV-2, obtained from 3 patients in Hong Kong that had a 5'-untranslated region segment at the end of the ORF6 gene that in its new location translated into an ORF6 protein with a predicted modified carboxyl terminus. While comparing the amino acid sequences of translated ORF8 genes in the GenBank database, we found a subsegment of the same 5'-UTR-derived amino acid sequence modifying the distal end of ORF8 of an isolate from the United States and decided to carry out a systematic search. METHODS: Using the nucleotide and in the case of SARS-CoV-2 also the translated amino acid sequence in three reading frames of the genomic termini of coronaviruses as query sequences, we searched for 5'-UTR sequences in regions other than the 5'-UTR in SARS-CoV-2 and reference strains of alpha-, beta-, gamma-, and delta-coronaviruses. RESULTS: We here report numerous genomic insertions of 5'-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses, but not delta- or gammacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions would change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain inserted 5'-UTR sequences. CONCLUSION: The intragenomic rearrangements involving 5'-untranslated region sequences described here, which in several cases affect highly conserved genes with a low propensity for recombination, may underlie the generation of variants homotypic with those of concern or interest and with potentially differing pathogenic profiles. Intragenomic rearrangements thus add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

Optimization of 5'UTR to evade SARS-CoV-2 Nonstructural protein 1-directed inhibition of protein synthesis in cells.

Maximizing the expression level of therapeutic proteins in cells is the general goal for DNA/mRNA therapies. It is particularly challenging to achieve efficient protein expression in the cellular contexts with inhibited translation machineries, such as in the presence of cellular Nonstructural protein 1 (Nsp1) of coronaviruses (CoVs) that has been reported to inhibit overall protein synthesis of host genes and exogenously delivered mRNAs/DNAs. In this study, we thoroughly examined the sequence and structure contexts of viral and non-viral 5'UTRs that determine the protein expression levels of exogenously delivered DNAs and mRNAs in cells expressing SARS-CoV-2 Nsp1. It was found that high 5'-proximal A/U content promotes an escape from Nsp1-directed inhibition of protein synthesis and results in selective protein expression. Furthermore, 5'-proximal Cs were found to significantly enhance the protein expression in an Nsp1-dependent manner, while Gs located at a specific window close to the 5'-end counteract such enhancement. The distinct protein expression levels resulted from different 5'UTRs were found correlated to Nsp1-induced mRNA degradations. These findings ultimately enabled rational designs for optimized 5'UTRs that lead to strong expression of exogenous proteins regardless of the translationally repressive Nsp1. On the other hand, we have also identified several 5'-proximal sequences derived from host genes that are capable of mediating the escapes. These results provided novel perspectives to the optimizations of 5'UTRs for DNA/mRNA therapies and/or vaccinations, as well as shedding light on the potential host escapees from Nsp1-directed translational shutoffs. KEY POINTS: • The 5'-proximal SL1 and 5a/b derived from SARS-CoV-2 genomic RNA promote exogenous protein synthesis in cells expressing Nsp1 comparing with non-specific 5'UTRs. • Specific 5'-proximal sequence contexts are the key determinants of the escapes from Nsp1-directed translational repression and thereby enhance protein expressions. • Systematic mutagenesis identified optimized 5'UTRs that strongly enhance protein expression and promote resistance to Nsp1-induced translational repression and RNA degradation.

eIF3d: A driver of noncanonical cap-dependent translation of specific mRNAs and a trigger of biological/pathological processes.

Eukaryotic initiation factor 3d (eIF3d), a known RNA-binding subunit of the eIF3 complex, is a 66 to 68-kDa protein with an RNA-binding motif and a cap-binding domain. Compared with other eIF3 subunits, eIF3d is relatively understudied. However, recent progress in studying eIF3d has revealed a number of intriguing findings on its role in maintaining eIF3 complex integrity, global protein synthesis, and in biological and pathological processes. It has also been reported that eIF3d has noncanonical functions in regulating translation of a subset of mRNAs by binding to 5'-UTRs or interacting with other proteins independent of the eIF3 complex and additional functions in regulating protein stability. The noncanonical regulation of mRNA translation or protein stability may contribute to the role of eIF3d in biological processes such as metabolic stress adaptation and in disease onset and progression including severe acute respiratory syndrome coronavirus 2 infection, tumorigenesis, and acquired immune deficiency syndrome. In this review, we critically evaluate the recent studies on these aspects of eIF3d and assess prospects in understanding the function of eIF3d in regulating protein synthesis and in biological and pathological processes.

The 5'UTR of HCoV-OC43 adopts a topologically constrained structure to intrinsically repress translation.

The emergence of SARS-CoV-2, which is responsible for the COVID-19 pandemic, has highlighted the need for rapid characterization of viral mechanisms associated with cellular pathogenesis. Viral UTRs represent conserved genomic elements that contribute to such mechanisms. Structural details of most CoV UTRs are not available, however. Experimental approaches are needed to allow for the facile generation of high-quality viral RNA tertiary structural models, which can facilitate comparative mechanistic efforts. By integrating experimental and computational techniques, we herein report the efficient characterization of conserved RNA structures within the 5'UTR of the HCoV-OC43 genome, a lab-tractable model coronavirus. We provide evidence that the 5'UTR folds into a structure with well-defined stem-loops (SLs) as determined by chemical probing and direct detection of hydrogen bonds by NMR. We combine experimental base-pair restraints with global structural information from SAXS to generate a 3D model that reveals that SL1-4 adopts a topologically constrained structure wherein SLs 3 and 4 coaxially stack. Coaxial stacking is mediated by short linker nucleotides and allows SLs 1 to 2 to sample different cojoint orientations by pivoting about the SL3,4 helical axis. To evaluate the functional relevance of the SL3,4 coaxial helix, we engineered luciferase reporter constructs harboring the HCoV-OC43 5'UTR with mutations designed to abrogate coaxial stacking. Our results reveal that the SL3,4 helix intrinsically represses translation efficiency since the destabilizing mutations correlate with increased luciferase expression relative to wildtype without affecting reporter mRNA levels, thus highlighting how the 5'UTR structure contributes to the viral mechanism.

Time Series Analysis of SARS-CoV-2 Genomes and Correlations among Highly Prevalent Mutations.

The efforts of the scientific community to tame the recent pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) seem to have been diluted by the emergence of new viral strains. Therefore, it is imperative to understand the effect of mutations on viral evolution. We performed a time series analysis on 59,541 SARS-CoV-2 genomic sequences from around the world to gain insights into the kinetics of the mutations arising in the viral genomes. These 59,541 genomes were grouped according to month (January 2020 to March 2021) based on the collection date. Meta-analysis of these data led us to identify significant mutations in viral genomes. Pearson correlation of these mutations led us to the identification of 16 comutations. Among these comutations, some of the individual mutations have been shown to contribute to viral replication and fitness, suggesting a possible role of other unexplored mutations in viral evolution. We observed that the mutations 241C>T in the 5' untranslated region (UTR), 3037C>T in nsp3, 14408C>T in the RNA-dependent RNA polymerase (RdRp), and 23403A>G in spike are correlated with each other and were grouped in a single cluster by hierarchical clustering. These mutations have replaced the wild-type nucleotides in SARS-CoV-2 sequences. Additionally, we employed a suite of computational tools to investigate the effects of T85I (1059C>T), P323L (14408C>T), and Q57H (25563G>T) mutations in nsp2, RdRp, and the ORF3a protein of SARS-CoV-2, respectively. We observed that the mutations T85I and Q57H tend to be deleterious and destabilize the respective wild-type protein, whereas P323L in RdRp tends to be neutral and has a stabilizing effect. IMPORTANCE We performed a meta-analysis on SARS-CoV-2 genomes categorized by collection month and identified several significant mutations. Pearson correlation analysis of these significant mutations identified 16 comutations having absolute correlation coefficients of >0.4 and a frequency of >30% in the genomes used in this study. The correlation results were further validated by another statistical tool called hierarchical clustering, where mutations were grouped in clusters on the basis of their similarity. We identified several positive and negative correlations among comutations in SARS-CoV-2 isolates from around the world which might contribute to viral pathogenesis. The negative correlations among some of the mutations in SARS-CoV-2 identified in this study warrant further investigations. Further analysis of mutations such as T85I in nsp2 and Q57H in ORF3a protein revealed that these mutations tend to destabilize the protein relative to the wild type, whereas P323L in RdRp is neutral and has a stabilizing effect. Thus, we have identified several comutations which can be further characterized to gain insights into SARS-CoV-2 evolution.

Translation of SARS-CoV-2 gRNA Is Extremely Efficient and Competitive despite a High Degree of Secondary Structures and the Presence of an uORF.

The SARS-CoV-2 infection generates up to nine different sub-genomic mRNAs (sgRNAs), in addition to the genomic RNA (gRNA). The 5'UTR of each viral mRNA shares the first 75 nucleotides (nt.) at their 5'end, called the leader, but differentiates by a variable sequence (0 to 190 nt. long) that follows the leader. As a result, each viral mRNA has its own specific 5'UTR in term of length, RNA structure, uORF and Kozak context; each one of these characteristics could affect mRNA expression. In this study, we have measured and compared translational efficiency of each of the ten viral transcripts. Our data show that most of them are very efficiently translated in all translational systems tested. Surprisingly, the gRNA 5'UTR, which is the longest and the most structured, was also the most efficient to initiate translation. This property is conserved in the 5'UTR of SARS-CoV-1 but not in MERS-CoV strain, mainly due to the regulation imposed by the uORF. Interestingly, the translation initiation mechanism on the SARS-CoV-2 gRNA 5'UTR requires the cap structure and the components of the eIF4F complex but showed no dependence in the presence of the poly(A) tail in vitro. Our data strongly suggest that translation initiation on SARS-CoV-2 mRNAs occurs via an unusual cap-dependent mechanism.

Cap-independent translation and a precisely located RNA sequence enable SARS-CoV-2 to control host translation and escape anti-viral response.

Translation of SARS-CoV-2-encoded mRNAs by the host ribosomes is essential for its propagation. Following infection, the early expressed viral protein NSP1 binds the ribosome, represses translation, and induces mRNA degradation, while the host elicits an anti-viral response. The mechanisms enabling viral mRNAs to escape this multifaceted repression remain obscure. Here we show that expression of NSP1 leads to destabilization of multi-exon cellular mRNAs, while intron-less transcripts, such as viral mRNAs and anti-viral interferon genes, remain relatively stable. We identified a conserved and precisely located cap-proximal RNA element devoid of guanosines that confers resistance to NSP1-mediated translation inhibition. Importantly, the primary sequence rather than the secondary structure is critical for protection. We further show that the genomic 5'UTR of SARS-CoV-2 drives cap-independent translation and promotes expression of NSP1 in an eIF4E-independent and Torin1-resistant manner. Upon expression, NSP1 further enhances cap-independent translation. However, the sub-genomic 5'UTRs are highly sensitive to eIF4E availability, rendering viral propagation partially sensitive to Torin1. We conclude that the combined NSP1-mediated degradation of spliced mRNAs and translation inhibition of single-exon genes, along with the unique features present in the viral 5'UTRs, ensure robust expression of viral mRNAs. These features can be exploited as potential therapeutic targets.

Specific Recognition of the 5'-Untranslated Region of West Nile Virus Genome by Human Innate Immune System.

In the last few years, the sudden outbreak of COVID-19 caused by SARS-CoV-2 proved the crucial importance of understanding how emerging viruses work and proliferate, in order to avoid the repetition of such a dramatic sanitary situation with unprecedented social and economic costs. West Nile Virus is a mosquito-borne pathogen that can spread to humans and induce severe neurological problems. This RNA virus caused recent remarkable outbreaks, notably in Europe, highlighting the need to investigate the molecular mechanisms of its infection process in order to design and propose efficient antivirals. Here, we resort to all-atom Molecular Dynamics simulations to characterize the structure of the 5'-untranslated region of the West Nile Virus genome and its specific recognition by the human innate immune system via oligoadenylate synthetase. Our simulations allowed us to map the interaction network between the viral RNA and the host protein, which drives its specific recognition and triggers the host immune response. These results may provide fundamental knowledge that can assist further antivirals' design, including therapeutic RNA strategies.

Evidence for selection on SARS-CoV-2 RNA translation revealed by the evolutionary dynamics of mutations in UTRs and CDSs.

RNA translation is the rate-limiting step when cells synthesize proteins. Elevating translation efficiency (TE) is intuitively beneficial. Particularly, when viruses invade host cells, how to compete with endogenous RNAs for efficient translation is a major issue to be resolved. We collected millions of worldwide SARS-CoV-2 sequences during the past year and traced the dynamics of allele frequency of every mutation. We defined adaptive and deleterious mutations according to the rise and fall of their frequencies along time. For 5'UTR and synonymous mutations in SARS-CoV-2, the selection on TE is evident near start codons. Adaptive mutations generally decrease GC content while deleterious mutations increase GC content. This trend fades away with increasing distance to start codons. Mutations decreasing GC content near start codons would unravel the complex RNA structure and facilitate translation initiation, which are beneficial to SARS-CoV-2, and vice versa. During this evolutionary arms race between human and virus, SARS-CoV-2 tries to improve its cis elements to compete with host RNAs for rapid translation.

1H, 13C and 15N chemical shift assignment of the stem-loops 5b + c from the 5'-UTR of SARS-CoV-2.

The ongoing pandemic of the respiratory disease COVID-19 is caused by the SARS-CoV-2 (SCoV2) virus. SCoV2 is a member of the Betacoronavirus genus. The 30 kb positive sense, single stranded RNA genome of SCoV2 features 5'- and 3'-genomic ends that are highly conserved among Betacoronaviruses. These genomic ends contain structured cis-acting RNA elements, which are involved in the regulation of viral replication and translation. Structural information about these potential antiviral drug targets supports the development of novel classes of therapeutics against COVID-19. The highly conserved branched stem-loop 5 (SL5) found within the 5'-untranslated region (5'-UTR) consists of a basal stem and three stem-loops, namely SL5a, SL5b and SL5c. Both, SL5a and SL5b feature a 5'-UUUCGU-3' hexaloop that is also found among Alphacoronaviruses. Here, we report the extensive 1H, 13C and 15N resonance assignment of the 37 nucleotides (nts) long sequence spanning SL5b and SL5c (SL5b + c), as basis for further in-depth structural studies by solution NMR spectroscopy.

Rocaglates as Antivirals: Comparing the Effects on Viral Resistance, Anti-Coronaviral Activity, RNA-Clamping on eIF4A and Immune Cell Toxicity.

Rocaglates are potent broad-spectrum antiviral compounds with a promising safety profile. They inhibit viral protein synthesis for different RNA viruses by clamping the 5'-UTRs of mRNAs onto the surface of the RNA helicase eIF4A. Apart from the natural rocaglate silvestrol, synthetic rocaglates like zotatifin or CR-1-31-B have been developed. Here, we compared the effects of rocaglates on viral 5'-UTR-mediated reporter gene expression and binding to an eIF4A-polypurine complex. Furthermore, we analyzed the cytotoxicity of rocaglates on several human immune cells and compared their antiviral activities in coronavirus-infected cells. Finally, the potential for developing viral resistance was evaluated by passaging human coronavirus 229E (HCoV-229E) in the presence of increasing concentrations of rocaglates in MRC-5 cells. Importantly, no decrease in rocaglate-sensitivity was observed, suggesting that virus escape mutants are unlikely to emerge if the host factor eIF4A is targeted. In summary, all three rocaglates are promising antivirals with differences in cytotoxicity against human immune cells, RNA-clamping efficiency, and antiviral activity. In detail, zotatifin showed reduced RNA-clamping efficiency and antiviral activity compared to silvestrol and CR-1-31-B, but was less cytotoxic for immune cells. Our results underline the potential of rocaglates as broad-spectrum antivirals with no indications for the emergence of escape mutations in HCoV-229E.

Correlated sequence signatures are present within the genomic 5'UTR RNA and NSP1 protein in coronaviruses.

The 5'UTR part of coronavirus genomes plays key roles in the viral replication cycle and translation of viral mRNAs. The first 75-80 nt, also called the leader sequence, are identical for genomic mRNA and subgenomic mRNAs. Recently, it was shown that cooperative actions of a 5'UTR segment and the nonstructural protein NSP1 are essential for both the inhibition of host mRNAs and for specific translation of viral mRNAs. Here, sequence analyses of both the 5'UTR RNA segment and the NSP1 protein have been done for several coronaviruses, with special attention to the betacoronaviruses. The conclusions are: (i) precise specific molecular signatures can be found in both the RNA and the NSP1 protein; (ii) both types of signatures correlate between each other. Indeed, definite sequence motifs in the RNA correlate with sequence motifs in the protein, indicating a coevolution between the 5'UTR and NSP1 in betacoronaviruses. Experimental mutational data on 5'UTR and NSP1 from SARS-CoV-2 using cell-free translation extracts support these conclusions and show that some conserved key residues in the amino-terminal half of the NSP1 protein are essential for evasion to the inhibitory effect of NSP1 on translation.

Targeting stem-loop 1 of the SARS-CoV-2 5' UTR to suppress viral translation and Nsp1 evasion.

SARS-CoV-2 is a highly pathogenic virus that evades antiviral immunity by interfering with host protein synthesis, mRNA stability, and protein trafficking. The SARS-CoV-2 nonstructural protein 1 (Nsp1) uses its C-terminal domain to block the messenger RNA (mRNA) entry channel of the 40S ribosome to inhibit host protein synthesis. However, how SARS-CoV-2 circumvents Nsp1-mediated suppression for viral protein synthesis and if the mechanism can be targeted therapeutically remain unclear. Here, we show that N- and C-terminal domains of Nsp1 coordinate to drive a tuned ratio of viral to host translation, likely to maintain a certain level of host fitness while maximizing replication. We reveal that the stem-loop 1 (SL1) region of the SARS-CoV-2 5' untranslated region (5' UTR) is necessary and sufficient to evade Nsp1-mediated translational suppression. Targeting SL1 with locked nucleic acid antisense oligonucleotides inhibits viral translation and makes SARS-CoV-2 5' UTR vulnerable to Nsp1 suppression, hindering viral replication in vitro at a nanomolar concentration, as well as providing protection against SARS-CoV-2-induced lethality in transgenic mice expressing human ACE2. Thus, SL1 allows Nsp1 to switch infected cells from host to SARS-CoV-2 translation, presenting a therapeutic target against COVID-19 that is conserved among immune-evasive variants. This unique strategy of unleashing a virus' own virulence mechanism against itself could force a critical trade-off between drug resistance and pathogenicity.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2.

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

Structure of the 5' untranslated region in SARS-CoV-2 genome and its specific recognition by innate immune system via the human oligoadenylate synthase 1.

2'-5'-Oligoadenylate synthetase 1 (OAS1) is one of the key enzymes driving the innate immune system response to SARS-CoV-2 infection whose activity has been related to COVID-19 severity. OAS1 is a sensor of endogenous RNA that triggers the 2'-5'-oligoadenylate/RNase L pathway. Upon SARS-CoV-2 infection, OAS1 is responsible for the recognition of viral RNA and has been shown to possess a particularly high sensitivity for the 5'-untranslated (5'-UTR) RNA region, which is organized in a double-strand stem loop motif (SL1). Here we report the structure of the SL1/OAS1 complex also rationalizing the high affinity for OAS1.

Putative Secondary Structure at 5'UTR as a Potential Antiviral Target against SARS-CoV-2.

SARS-CoV-2 is an enveloped positive-sense single-stranded RNA coronavirus that causes COVID-19, of which the current outbreak has resulted in a high number of cases and fatalities throughout the world, even vaccine doses are being administered. The aim of this work was to scan the SARS-CoV-2 genome in search for therapeutic targets. We found a sequence in the 5'UTR (NC\_045512:74-130), consisting of a typical heptamer next to a structured region that may cause ribosomal frameshifting. The potential biological value of this region is relevant through its low similarity with other viruses, including coronaviruses related to SARS-CoV, and its high sequence conservation within multiple SARS-CoV-2 isolates. We have predicted the secondary structure of the region by means of different bioinformatic tools. We have suggested a most probable secondary structure to proceed with a 3D reconstruction of the structured segment. Finally, we carried out virtual docking on the 3D structure to look for a binding site and then for drug ligands from a database of lead compounds. Several molecules that could be probably administered as oral drugs show promising binding affinity within the structured region, and so it could be possible interfere its potential regulatory role.

Inhibition of SARS-CoV-2 coronavirus proliferation by designer antisense-circRNAs.

Circular RNAs (circRNAs) are noncoding RNAs that exist in all eukaryotes investigated and are derived from back-splicing of certain pre-mRNA exons. Here, we report the application of artificial circRNAs designed to act as antisense-RNAs. We systematically tested a series of antisense-circRNAs targeted to the SARS-CoV-2 genome RNA, in particular its structurally conserved 5'-untranslated region. Functional assays with both reporter transfections as well as with SARS-CoV-2 infections revealed that specific segments of the SARS-CoV-2 5'-untranslated region can be efficiently accessed by specific antisense-circRNAs, resulting in up to 90% reduction of virus proliferation in cell culture, and with a durability of at least 48 h. Presenting the antisense sequence within a circRNA clearly proved more efficient than in the corresponding linear configuration and is superior to modified antisense oligonucleotides. The activity of the antisense-circRNA is surprisingly robust towards point mutations in the target sequence. This strategy opens up novel applications for designer circRNAs and promising therapeutic strategies in molecular medicine.

Cross-sectional genomic perspective of epidemic waves of SARS-CoV-2: A pan India study.

BACKGROUND: COVID-19 has posed unforeseen circumstances and throttled major economies worldwide. India has witnessed two waves affecting around 31 million people representing 16% of the cases globally. To date, the epidemic waves have not been comprehensively investigated to understand pandemic progress in India. OBJECTIVE: Here, we aim for pan Indian cross-sectional evolutionary analysis since inception of SARS-CoV-2. METHODS: High quality genomes, along with their collection date till 26th July 2021, were downloaded. Whole genome-based phylogeny was obtained. Further, the mutational analysis was performed using SARS-CoV-2 first reported from Wuhan (NC_045512.2) as reference. RESULTS: Based on reported cases and mutation rates, we could divide the Indian epidemic into seven phases. The average mutation rate for the pre-first wave was <11, which elevated to 17 in the first wave and doubled in the second wave (∼34). In accordance with mutation rate, VOCs and VOIs started appearing in the first wave (1.5%), which dominated the second (∼96%) and post-second wave (100%). Nation-wide mutational analysis depicted >0.5 million mutation events with four major mutations in >19,300 genomes, including two mutations in coding (spike (D614G), and NSP 12b (P314L) of rdrp), one silent mutation (NSP3 F106F) and one extragenic mutation (5' UTR 241). CONCLUSION: Whole genome-based phylogeny could demarcate post-first wave isolates from previous ones by point of diversification leading to incidences of VOCs and VOIs in India. Such analysis is crucial in the timely management of pandemic.

Insights on the SARS-CoV-2 genome variability: the lesson learned in Brazil and its impacts on the future of pandemics.

Since the beginning of the SARS-CoV-2 spread in Brazil, few studies have been published analysing the variability of viral genome. Herein, we described the dynamic of SARS-CoV-2 strains circulating in Brazil from May to September 2020, to better understand viral changes that may affect the ongoing pandemic. Our data demonstrate that some of the mutations identified are currently observed in variants of interest and variants of concern, and emphasize the importance of studying previous periods in order to comprehend the emergence of new variants. From 720 SARS-CoV-2 genome sequences, we found few sites under positive selection pressure, such as the D614G (98.5 %) in the spike, that has replaced the old variant; the V1167F in the spike (41 %), identified in the P.2 variant that emerged from Brazil during the period of analysis; and I292T (39 %) in the N protein. There were a few alterations in the UTRs, which was expected, however, our data suggest that the emergence of new variants was not influenced by mutations in UTR regions, since it maintained its conformational structure in most analysed sequences. In phylogenetic analysis, the spread of SARS-CoV-2 from the large urban centres to the countryside during these months could be explained by the flexibilization of social isolation measures and also could be associated with possible new waves of infection. These results allow a better understanding of SARS-CoV-2 strains that have circulated in Brazil, and thus, with relevant infomation, provide the potential viral changes that may have affected and/or contributed to the current and future scenario of the COVID-19 pandemic.